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In an experiment designed to search for and study very rare two-body 
decay modes of the K®, we have observed four examples of the decay K\ — > 
e + e~, where the expected background is 0.17 ±0.10 events. This observation 
translates into a branching fraction of 8.7~^\'\ x 10 -12 , consistent with recent 
theoretical predictions. This result represents by far the smallest branching 
fraction yet measured in particle physics. 
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We report the first observation of the decay K\ — > e + e~, which is strongly suppressed 
by the GIM mechanism [|l| and the helicity structure of the V — A interaction. Thus, it 
could be sensitive to a new physics process not subject to these effects. The rate of the 
analogous decay K° L — > /x 4 pT is dominated by the absorptive contribution which arises from 
the real two-photon intermediate state, K\ — > 77 — > The absorptive contribution |2j 

suggests a K\ — > e + e~ branching fraction of 3xl0 -12 . However, recent predictions || in the 
framework of chiral perturbation theory indicate that the two-photon dispersive contribution 
is larger than the absorptive part in the K® — *> e + e~ decay. The predicted branching fraction 
is approximately 9 x 10~ 12 . The observation reported here is consistent with these predictions 
and consequently rules out a significant non-Standard Model contribution to this decay. 

Data were recorded in 1995 and 1996 by Experiment 871 at the Alternating Gradient 
Synchrotron (AGS) at Brookhaven National Laboratory (BNL). The experiment built upon 
experience with BNL E791, which set the best upper limit on K\ — *> e + e~ [|J. An improved 
search was made possible by the large increase in AGS intensity with the new Booster and 
by apparatus modifications to take advantage of higher flux while improving background 
rejection. Major parts of the spectrometer were more finely segmented to accommodate 
higher rates. Other important changes were: (a) redundancy was increased in the most 
critical areas by making three position measurements in all tracking stations; (b) magnetic 
fields were set to provide tracks nearly parallel to the beam for two-body decays at the 
Jacobian peak; and (c) a beam stop || was installed in the first spectrometer magnet to 
absorb the neutral beam. Taken together, these modifications made possible an improvement 
of roughly an order of magnitude over the earlier experiment. 

The neutral beam was produced by 24 GeV protons on a 1.4 interaction length Pt target 
at 3.75° with respect to the collimation channel. Typical proton intensity was 1.5 x 10 13 per 
slow spill (1.2-1.6 s duration). This provided about 2x 10 8 K\ per spill (2 < px < 16 GeV/c) 
at the entrance of the decay volume, of which 7.5% decayed in the 11m fiducial length. The 
experiment is shown in Fig. 0. For brevity, only the detectors most relevant to this analysis 
will be discussed. These are tracking, electron identification, and trigger hodoscopes. 

The tracking section consisted of six chamber stations and two consecutive dipole mag- 
nets. The magnets had opposite polarities and provided net transverse momentum kicks of 
418 MeV/c and 216 MeV/c. The topology of a K\ — > e + e~ decay is illustrated in Fig. [l|. 
The upstream four tracking stations, where the highest rates occurred, were constructed 
from 5 mm diameter straw tubes || and were operated with a fast gas mixture (50-50 
CF4/C2H 6 ). The downstream-most two stations were in a region of lower rates and con- 
sisted of drift chambers with 1 cm sense wire spacing operated with a 50-50 Ar/C2Hg gas 
mixture. Each tracking station provided three measurements of position in the horizontal 
(magnet-bend) x-plane and two measurements in the vertical y-plane. Since the resolution 
of the system was limited mainly by multiple Coulomb scattering, the total material was 
kept to a minimum; including helium bags between the chambers, it amounted to 1.5 x 10 -2 
radiation lengths. The single wire resolutions were 160 and 120 /im, and the average effi- 
ciencies (per wire) were 96% and 98%, for the straw and drift chambers, respectively. The 
measured mass resolutions for K\ — > -n + n~ and K\ — > yU + /i~ decay modes were 1.11 and 
1.28 MeV/c 2 , respectively, compared with the 1.07 and 1.22 MeV/c 2 predicted by Monte 
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Carlo. The K\ — > e + e~ mass resolution predicted by Monte Carlo is 1.39 MeV/c 2 . 

Redundant electron identification was achieved by an atmospheric hydrogen threshold 
Cherenkov counter (CER) and a lead glass array (PBG). The CER had 4-by-4 arrays of 
mirror-phototube pairs on each side of the beam. The average photoelectron yield was 5.5 
for electrons. The PBG consisted of 216 blocks arranged in two layers, with 3.5 radiation 
length converter blocks in front of 10.5 radiation length absorbe r blocks. The measured 
PBG energy resolution for electrons was a/E = 0.015 + 0.062/y£^(GeV). The CER and 
PBG performance is summarized in Table |. 

Two scintillator hodoscopes, one upstream of the CER and one downstream (separated 
by 2.9 m), were used for triggering. Both hodoscopes had 3.2 cm wide x-measuring slats 
with phototubes on both ends. The downstream hodoscope had, in addition, 3.0 cm wide 
y-measuring slats with a phototube on one end. The slats overlapped their nearest neighbors 
by 3 mm to avoid inefficiency due to cracks. 

Online data selection involved hardware and software triggers. A Level (L0) trigger 
based on the pattern of trigger hodoscope hits required two in-time tracks to be parallel 
within about 30 mrad to, and on opposite sides of, the beam axis. The L0 trigger along with 
signals from particle identification detectors formed Level 1 (LI) triggers for various di-lepton 
modes. The LI trigger for e + e~ required CER hits to be in time and spatially correlated 
with the trigger hodoscope hits. PBG signals were not required. For events passing a 
LI trigger, a software algorithm, run on a farm of eight RISC processors, constituted a 
Level 3 (L3) trigger. The algorithm reconstructed tracks using information from the trigger 
hodoscopes and all tracking chambers. For e + e~ events, it required tracks on each side of 
the spectrometer which formed a vertex in the decay volume. Also, L0 triggers prescaled by 
a factor of 1000 were recorded for analysis. This minimum bias sample was used to study 
detector performance and provided data for normalization. 

The spectrometer had a geometrical acceptance of 1.89% for — > e + e~ decays with 
9.75 < z < 20.75 m (see Fig. 1) and kaon momenta between 2 and 16 GeV/c. The trig- 
ger requirement of parallelism reduced the acceptance to 1.57%. Analysis criteria further 
reduced the acceptance to 1.23%. 

Potential sources of background include misidentified K\ — > 7r ± e =F z/ decays, accidental 
spatial and temporal coincidences of e + and e~ from two K\ — ► n^e^v decays, partially 
measured K\ — > e + e~7, K\ — > e + e~e + e - , or K\ — > e + e~77 decays, and K\ — > 77 decays 
with asymmetric external conversion of both photons in the vacuum window or first straw 
chamber. The only sources which were not negligible after all analysis criteria were applied 
are K\ — > e + e~j J7[ and K\ — > e + e~e + e~ || decays. Both of these decays have a low 
probability of producing an e + e~ pair with invariant mass near the K\ mass. 

Analysis procedures were designed to ensure that selection criteria were not influenced by 
knowledge of events in or near the signal region. K\ — > e + e~ events should have measured 
e + e~ mass near the K\ mass (497.7 MeV/c 2 ) and measured transverse momentum (pr), 
defined with respect to the parent K\ direction of flight, near zero. Thus, events with 
490 < M ee < 505 MeV/c 2 and p\ < 100 (MeV/c) 2 were excluded from consideration until 
all selection criteria were finalized. 

As a first step in the analysis, events were required to have a good track (with signals 
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from at least two x-measuring wires and one ^-measuring wire in each chamber) on each side 
of the spectrometer. These tracks had to form a good vertex in the decay volume and project 
to in-time trigger counter hits consistent with the parallelism requirement. Two independent 
fitting algorithms, which used a full magnetic field map, subsequently determined kinematic 
quantities for each track. The fitters had different sensitivities to track finding errors and 
checked each other's results. Ultimately, all selection criteria were applied to the results of 
one fitter, with the exception that consistent results for M ee and p\ were required of both. 

Events had to satisfy minimal track and vertex quality criteria, be fully contained within 
the spectrometer, and have M ee > 475 MeV/c 2 and p\ < 900 (MeV/c) 2 (excluding the 
region 490 < M ee < 505 MeV/c 2 and p% < 100 (MeV/c) 2 ). In addition, events had to have 
energy deposition in the PBG characteristic of electrons and good CER pulse height and 
timing information. The remaining 833 events were studied to determine criteria to further 
exclude background with minimal acceptance loss. This optimization was based on Monte 
Carlo and studies of observed K® — > 7r + 7r~ events. These criteria included better track and 
vertex quality, tighter timing, a requirement that the momentum asymmetry between the 
tracks (\p e + — p e -\/(Pe+ +Pe-)) be less than 0.55 (to suppress K\ — > 7r ± e =F z/ decays), and a 
requirement that no additional complete tracks be found in the spectrometer (to suppress 
events with two K\ decays). Finally, to reduce K® — > e + e~e + e~ background, criteria were 
applied to reject events with short partial tracks in the two upstream-most straw chambers 
pointing to the reconstructed e + e~ vertex. The resulting sample comprised 44 events. 

After all analysis criteria were finalized, a signal region was defined (elliptical in p\ and 
M ee and corresponding to about 2.5 sigma in each). As a result of inner bremsstrahlung ||, 
23% of K\ — > e + e~ decays fall outside this ellipse. The size of the signal region was 
chosen to reduce the expected level of background to well below one event. Background 
estimates were based on Monte Carlo simulations and comparisons to the data. Samples of 
K\ — > e + e~e + e~ and K\ — > e + e~7 decays were generated [flQ | . These distributions were 
absolutely normalized by their measured branching fractions |L1] to obtain predictions of 
38 ± 8 events from K\ — > e + e~7 and 24 ± 11 events from K\ — > e + e~e + e~ in the region 
defined by 476 < M ee < 490 MeV/c 2 andp 2 . < 400 (MeV/c) 2 , where 43 events were observed. 
The uncertainties in the predictions are mainly due to the if£ — > e + e~7 form factor and 
the fact that the K\ — > e + e~e + e~ efficiency of the partial track cut is not well known. 
To remove these sources of systematic uncertainty from the estimate of the background in 
the signal region, we performed a fit to data to normalize the background distributions, 
which were then extrapolated into the signal region. The fit in p\ and M ee exploited the 
fact that in the range 476 < M ee < 490 MeV/c the p\ distributions for the background 
decays differ significantly. This procedure yielded estimates of background in the signal 
region of 0.09 ± 0.11 event from K\ — > e + e~e + e~, 0.08 ± 0.02 event from K\ — > e + e~7, 
and 0.17 ± 0.10 for the sum of both (to be compared with 0.37 ± 0.14 for the absolutely 
normalized prediction) . 

When the full analysis was performed on the excluded region, four K\ —>■ e + e~ candidates 
were found inside the signal ellipse, as shown in Fig. |2|. These candidates have been carefully 
scrutinized and exhibit no anomalous features. The probability of observing 4 background 
events in the signal region when 0.2 are expected is 6 x 10 -5 . A check on the reconstructed 
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mass of K\ — > 7i + ii~ events collected within minutes of the K\ — > e + e~ candidates rules 
out transient shifts in the mass scale. 

To extract the best estimate of the number of K\ — > e + e~ events, a maximum likelihood 
fit was performed in the region defined by 476 < M ee < 510 MeV/c 2 andp 2 - < 400 (MeV/c) 2 . 
In addition to the 50 events inside this region, the inputs to the fit were Monte Carlo 
distributions in M ee and p\ for K® — > e + e~e + e~ , K° — > e + e~7, and lf£ — > e + e~ decays. 
The fit estimated the number of events in each distribution, subject to the constraint that 
the sum equal the observed number of events. Fig. [3] shows the fit results for the three 
distributions versus M ee . The result for the number of K\ — > e + e~ events is 4.20l 2 ;g4 i n 
the full region of the fit. For comparison, the inset in Fig. |3| shows the M ee distribution 
of data along with the absolutely normalized Monte Carlo predictions for K\ — > e + e~e + e~ 
and K\ — > e + e _ 7. 

The K® — > e + e~ branching fraction B ee was determined from the formula 

-Dee = ^ee ~JtTt 4 7ET 7L3 TpIDJ™- 
ju, tit ^ee c ee c ee c ee 

N ee is the number of K\ — > e + e~ events determined from the likelihood fit. B nir is the K\ — > 
7r + 7r _ branching fraction [pT[| . R is the prescale factor for the minimum bias sample used for 



normalization. It was the product of the hardware prescale (1000) and an additional factor of 
20 imposed in software. N nn is the number of K\ — > tx + h~ events in this reduced sample (and 
includes a small correction for the effect of Kg interference). The ratio A n7T /A ee , calculated 
via Monte Carlo, corrects for mode-dependent acceptance differences due to the detector 
geometry, trigger, and analysis criteria, is the LI ee trigger efficiency with respect to the 
L0 trigger (which all events had to satisfy), eff is the L3 ee trigger efficiency and was derived 
from the measured L3 efficiency for K\ — > 7r + 7r~ by making a small correction to account 
for the different kinematics of K\ — ► e + e _ . e^J D is the particle identification efficiency and 
was determined by weighting measured detector efficiencies by the appropriate distributions 
for K° L — ► e + e~ decays. f nn accounts for the loss of K\ — > 7r + 7r~ events due to hadronic 
interactions in the spectrometer and was based on a GEANT calculation. The calculation 
was checked against special data taken with a single-arm trigger, which provided events in 
which a pion was lost due to hadronic interaction in the non-trigger arm. The values of these 
parameters are given in Table |TJ. The calculated branching fraction B ee is 8.7±|:I x 10- 12 . 

In summary, we have made the first observation of the decay K\ — ► e + e~. The measured 
branching fraction is consistent with recent theoretical predictions and is the smallest 
measured in particle physics to date. 
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TABLE I. Performance of particle identification detectors as measured from a clearly identified 
sample of K e ^ and events. The Cherenkov it and \x rejections are calculated for particles with 
momenta below their Cherenkov thresholds, 8.3 GeV/c and 6.3 GeV/c, respectively. 







Cherenkov 


Lead Glass 


e 

7T 

n 


efficiency 

misidentification 

misidentification 


0.977± 0.001 
0.0019 ± 0.0002 
0.0024 ± 0.0002 


0.987 ± 0.004 
0.0093 ± 0.0004 
0.0018 ± 0.0002 




TABLE II. Factors entering 


into the calculation of the — > e + e branching fraction. 




Variable 




Value 




N ee 

BirTr [11 

R 

N 

LI 

*-ee 

,L3 
•-ee 
PID 
t ee 

firir 




4.20l?j! 
(2.067 ± 0.035) x I0~ 3 
20000 
83531 ± 381 
1.478 ±0.011 
0.977 ±0.011 
0.929 ± 0.003 
0.929 ±0.011 
0.954 ± 0.004 
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FIG. 1. Plan view of the E871 apparatus. The origin of the z-axis is at the target. 
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FIG. 2. p?j, versus M ee for K\ e + e candidates. The dashed line shows the exclusion region. 
The solid curve bounds the signal region. 
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FIG. 3. Results of the maximum likelihood fit, showing the distributions of — > e + e~, 
K 1 ^ — > e + e~e + e _ , and if£ — > e + e~7 along with the data versus M ee . The inset shows the abso- 
lute Monte Carlo prediction for the background distributions, which have systematic uncertainties 
described in the text. 
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